A study of the sensitivity of shape measurements to the input parameters of weak lensing image simulations

Improvements in the accuracy of shape measurements are essential to exploit the statistical power of planned imaging surveys that aim to constrain cosmological parameters using weak lensing by large-scale structure. Although a range of tests can be performed using the measurements, the performance of the algorithm can only be quantified using simulated images. This yields, however, only meaningful results if the simulated images resemble the real observations sufficiently well. In this paper we explore the sensitivity of the multiplicative bias to the input parameters of Euclid-like image simulations.We find that algorithms will need to account for the local density of sources. In particular the impact of galaxies below the detection limit warrants further study, because magnification changes their number density, resulting in correlations between the lensing signal and multiplicative bias. Although achieving sub-percent accuracy will require further study, we estimate that sufficient archival Hubble Space Telescope data are available to create realistic populations of galaxies.Comment: 18 pages, accepted for publications in MNRA

Evolution of the luminosity-to-halo mass relation of LRGs from a combined SDSS-DR10+RCS2 analysis

We study the evolution of the luminosity-to-halo mass relation of Luminous Red Galaxies (LRGs). We select a sample of 52 000 LOWZ and CMASS LRGs from the Baryon Oscillation Spectroscopic Survey (BOSS) SDSS-DR10 in the ~450 deg^2 that overlaps with imaging data from the second Red-sequence Cluster Survey (RCS2), group them into bins of absolute magnitude and redshift and measure their weak lensing signals. The source redshift distribution has a median of 0.7, which allows us to study the lensing signal as a function of lens redshift. We interpret the lensing signal using a halo model, from which we obtain the halo masses as well as the normalisations of the mass-concentration relations. We find that the concentration of haloes that host LRGs is consistent with dark matter only simulations once we allow for miscentering or satellites in the modelling. The slope of the luminosity-to-halo mass relation has a typical value of 1.4 and does not change with redshift, but we do find evidence for a change in amplitude: the average halo mass of LOWZ galaxies increases by 25_{-14}^{+16} % between z=0.36 and 0.22 to an average value of 6.43+/-0.52 x 10^13 h70^-1 Msun. If we extend the redshift range using the CMASS galaxies and assume that they are the progenitors of the LOWZ sample, we find that the average mass of LRGs increases by 80^{+39}_{-28} % between z=0.6 and 0.2Comment: 20 pages, 11 figures, accepted for publication in A&

The Canadian Cluster Comparison Project: detailed study of systematics and updated weak lensing masses

Masses of clusters of galaxies from weak gravitational lensing analyses of ever larger samples are increasingly used as the reference to which baryonic scaling relations are compared. In this paper we revisit the analysis of a sample of 50 clusters studied as part of the Canadian Cluster Comparison Project. We examine the key sources of systematic error in cluster masses. We quantify the robustness of our shape measurements and calibrate our algorithm empirically using extensive image simulations. The source redshift distribution is revised using the latest state-of-the-art photometric redshift catalogs that include new deep near-infrared observations. Nonetheless we find that the uncertainty in the determination of photometric redshifts is the largest source of systematic error for our mass estimates. We use our updated masses to determine b, the bias in the hydrostatic mass, for the clusters detected by Planck. Our results suggest 1-b=0.76+-0.05(stat)}+-0.06(syst)}, which does not resolve the tension with the measurements from the primary cosmic microwave background.Comment: resubmitted to MNRAS after review by refere

Brightest Cluster Galaxies Trace Weak Lensing Mass Bias and Halo Triaxiality in The Three Hundred Project

Galaxy clusters have a triaxial matter distribution. The weak-lensing signal, an important part in cosmological studies, measures the projected mass of all matter along the line-of-sight, and therefore changes with the orientation of the cluster. Studies suggest that the shape of the brightest cluster galaxy (BCG) in the centre of the cluster traces the underlying halo shape, enabling a method to account for projection effects. We use 324 simulated clusters at four redshifts between 0.1 and 0.6 from `The Three Hundred Project' to quantify correlations between the orientation and shape of the BCG and the halo. We find that haloes and their embedded BCGs are aligned, with an average $\sim$20 degree angle between their major axes. The bias in weak lensing cluster mass estimates correlates with the orientation of both the halo and the BCG. Mimicking observations, we compute the projected shape of the BCG, as a measure of the BCG orientation, and find that it is most strongly correlated to the weak-lensing mass for relaxed clusters. We also test a 2-dimensional cluster relaxation proxy measured from BCG mass isocontours. The concentration of stellar mass in the projected BCG core compared to the total stellar mass provides an alternative proxy for the BCG orientation. We find that the concentration does not correlate to the weak-lensing mass bias, but does correlate with the true halo mass. These results indicate that the BCG shape and orientation for large samples of relaxed clusters can provide information to improve weak-lensing mass estimates.Comment: 15 pages, 9 figures, Figure 7 is key plot, Updated to match version accepted by MNRA

Towards emulating cosmic shear data:Revisiting the calibration of the shear measurements for the Kilo-Degree Survey

Exploiting the full statistical power of future cosmic shear surveys will necessitate improvements to the accuracy with which the gravitational lensing signal is measured. We present a framework for calibrating shear with image simulations that demonstrates the importance of including realistic correlations between galaxy morphology, size and more importantly, photometric redshifts. This realism is essential so that selection and shape measurement biases can be calibrated accurately for a tomographic cosmic shear analysis. We emulate Kilo-Degree Survey (KiDS) observations of the COSMOS field using morphological information from {\it Hubble} Space Telescope imaging, faithfully reproducing the measured galaxy properties from KiDS observations of the same field. We calibrate our shear measurements from lensfit, and find through a range of sensitivity tests that lensfit is robust and unbiased within the allowed 2 per cent tolerance of our study. Our results show that the calibration has to be performed by selecting the tomographic samples in the simulations, consistent with the actual cosmic shear analysis, because the joint distributions of galaxy properties are found to vary with redshift. Ignoring this redshift variation could result in misestimating the shear bias by an amount that exceeds the allowed tolerance. To improve the calibration for future cosmic shear analyses, it will be essential to also correctly account for the measurement of photometric redshifts, which requires simulating multi-band observations.Comment: 31 pages, 17 figures and 2 tables. Accepted for publication in A&A. Matches the published versio

KiDS-450: testing extensions to the standard cosmological model

We test extensions to the standard cosmological model with weak gravitational lensing tomography using 450 deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). In these extended cosmologies, which include massive neutrinos, nonzero curvature, evolving dark energy, modified gravity, and running of the scalar spectral index, we also examine the discordance between KiDS and cosmic microwave background measurements from Planck. The discordance between the two datasets is largely unaffected by a more conservative treatment of the lensing systematics and the removal of angular scales most sensitive to nonlinear physics. The only extended cosmology that simultaneously alleviates the discordance with Planck and is at least moderately favored by the data includes evolving dark energy with a time-dependent equation of state (in the form of the $w_0-w_a$ parameterization). In this model, the respective $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ constraints agree at the $1\sigma$ level, and there is `substantial concordance' between the KiDS and Planck datasets when accounting for the full parameter space. Moreover, the Planck constraint on the Hubble constant is wider than in LCDM and in agreement with the Riess et al. (2016) direct measurement of $H_0$. The dark energy model is moderately favored as compared to LCDM when combining the KiDS and Planck measurements, and remains moderately favored after including an informative prior on the Hubble constant. In both of these scenarios, marginalized constraints in the $w_0-w_a$ plane are discrepant with a cosmological constant at the $3\sigma$ level. Moreover, KiDS constrains the sum of neutrino masses to 4.0 eV (95% CL), finds no preference for time or scale dependent modifications to the metric potentials, and is consistent with flatness and no running of the spectral index. The analysis code is public at https://github.com/sjoudaki/kids450Comment: 22 pages, 16 figures, results unchanged, version accepted for publication by MNRA

The third data release of the Kilo-Degree Survey and associated data products

The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope. It aims to image 1500 square degrees in four filters (ugri). The core science driver is mapping the large-scale matter distribution in the Universe, using weak lensing shear and photometric redshift measurements. Further science cases include galaxy evolution, Milky Way structure, detection of high-redshift clusters, and finding rare sources such as strong lenses and quasars. Here we present the third public data release (DR3) and several associated data products, adding further area, homogenized photometric calibration, photometric redshifts and weak lensing shear measurements to the first two releases. A dedicated pipeline embedded in the Astro-WISE information system is used for the production of the main release. Modifications with respect to earlier releases are described in detail. Photometric redshifts have been derived using both Bayesian template fitting, and machine-learning techniques. For the weak lensing measurements, optimized procedures based on the THELI data reduction and lensfit shear measurement packages are used. In DR3 stacked ugri images, weight maps, masks, and source lists for 292 new survey tiles (~300 sq.deg) are made available. The multi-band catalogue, including homogenized photometry and photometric redshifts, covers the combined DR1, DR2 and DR3 footprint of 440 survey tiles (447 sq.deg). Limiting magnitudes are typically 24.3, 25.1, 24.9, 23.8 (5 sigma in a 2 arcsec aperture) in ugri, respectively, and the typical r-band PSF size is less than 0.7 arcsec. The photometric homogenization scheme ensures accurate colors and an absolute calibration stable to ~2% for gri and ~3% in u. Separately released are a weak lensing shear catalogue and photometric redshifts based on two different machine-learning techniques.Comment: small modifications; 27 pages, 12 figures, accepted for publication in Astronomy & Astrophysic

Shear nulling after PSF Gaussianisation: Moment-based weak lensing measurements with subpercent noise bias

Context. Current optical imaging surveys for cosmology cover large areas of sky. Exploiting the statistical power of these surveys for weak lensing measurements requires shape measurement methods with subpercent systematic errors. Aims. We introduce a new weak lensing shear measurement algorithm, shear nulling after PSF Gaussianisation (SNAPG), designed to avoid the noise biases that affect most other methods. Methods. SNAPG operates on images that have been convolved with a kernel that renders the point spread function (PSF) a circular Gaussian, and uses weighted second moments of the sources. The response of such second moments to a shear of the pre-seeing galaxy image can be predicted analytically, allowing us to construct a shear nulling scheme that finds the shear parameters for which the observed galaxies are consistent with an unsheared, isotropically oriented population of sources. The inverse of this nulling shear is then an estimate of the gravitational lensing shear. Results. We identify the uncertainty of the estimated centre of each galaxy as the source of noise bias, and incorporate an approximate estimate of the centroid covariance into the scheme. We test the method on extensive suites of simulated galaxies of increasing complexity, and find that it is capable of shear measurements with multiplicative bias below 0.5 percent